Static phase and dynamic scaling in a deposition model with an inactive species 
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We extend a previously proposed deposition model with two kinds of particles, considering the 
restricted solid-on-solid condition. The probability of incidence of particle C (A) is p (1 — p). 
Aggregation is possible if the top of the column of incidence has a nearest neighbor A and if the 
difference in the heights of neighboring columns does not exceed 1. For any value of p > 0, the 
deposit attains some static configuration, in which no deposition attempt is accepted. In 1 + 1 
dimensions, the interface width has a limiting value W s ~ P~ v , with rj — 3/2, which is confirmed 
by numerical simulations. The dynamic scaling relation W s = p~ v f (tp z ) is obtained in very large 
substrates, with z = r\. 

PACS numbers: 05.40.+j, 05.50.+q 
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I. INTRODUCTION 

Statistical growth models of surfaces and interfaces 
have attracted many attention in the last two decades 
motivated by technological applications of thin films and 
related nanostructuresoii In recent works, models with 
two types of particles were introduced, in order to repre- 
sent the effects .of-different chemical species in the depo- 
sition processesQ~li3. The competition of different growth 
mechanisms may lead to crossover of growth exponents 
and roughening transitions, as observed in many systems 
with a single specieal3iH3. 

A particularly interesting two-species model was pro- 
posed by Wang and CerdeiraO, which will be called AC 
model. In that model, particles A and C are released with 
probabilities 1—p and p, respectively, and aggregation is 
allowed only if the incident particle encounters a neigh- 
boring A at the sticking position (which may be defined 
by different rules). Thus, particles C represent impuri- 
ties that block the growth in their neighborhoods. For 
high p, the surface will be contaminated with this species 
and the growth process will fail. In previous works, the 
crossoyetTof growth exponents was studied in the growth 
regimeOu! 

In the present work, we will consider the restricted 
solid-on-solid (RSOS) version of the AC model. The 
RSOS model was introduced by Kim and Kosterlitz in 
1989 to describe the growth of thin films in which the 
heights' differences between neighboring columns do not 
exceed certain limiting value AH max . This condition 
prevents the formation of high local slopes in the film sur- 
faced, then it is interesting for the description of depo- 
sition processes in which diffusion and desorption mech- 
anisms (not explicitly included in the model) favor the 
formation of locally smooth surfaces. 

The RSOS version of the AC model is defined as fol- 
lows. At each deposition attempt, an incident particle, A 
or C, is chosen with the probabilities 1—p and p, respec- 
tively. This particle is released above a d-dimensional 



substrate in a randomly chosen column. The sticking 
position for the incident particle is the top of the se- 
lected column, but the aggregation is possible only if 
both the following conditions are satisfied: (a) the dif- 
ference in the heights of neighboring columns do not ex- 
ceed AH max = 1; (b) the sticking position has a nearest 
neighbor particle A. If one or both conditions are not 
satisfied, then the deposition attempt is rejected. Fig. 1 
illustrates the deposition rules. 
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FIG. 1. Examples of application of the deposition rules of 
the RSOS version of the AC model. Open squares represent 
particles A and filled squares represent particles C. In (a), 
the aggregation attempt is not accepted because there is no 
neighboring A at the top of the column. In (b), this neighbor 
is present (the dashed square indicates the sticking position). 
In (c) and (d), the aggregation attempt is not accepted be- 
cause it would violate the RSOS condition. 



Here we will study the model in d = 1. We will 
show that a dynamic transition occurs at p = because 
any finite flux of particles C will eventually suppress the 
growth process. Thus, at p > the model presents a 
static phase, i. c. the film attains a configuration that 
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cannot continue growing because no deposition attempt 
can be accepted. The interface width at saturation scales 
with p with an exponent r\ that can be exactly obtained. 
It is also shown that the dynamic exponent of the model 
is z = rj. The features of this static phase differ from 
the dynamic nature of the smooth phases of other mod- 
els with roughening transitions, such as those includ- 
ing competition between adsorption and desorption of 
adatomscij EJ. However, there are many important open 
questions in the field of roughening transitions, such as 
those concerning exponents' relationsli-3, then some re- 
sults presented here may also be helpful in that context. 

The rest of the paper is organized as follows. In Sec. II 
we present the simulations' results and discuss the tran- 
sition at p = 0. In Sec. Ill, we obtain a dynamic scaling 
relation for the interface width. In Sec. IV we present 
our conclusions. 



that the deposit is faceted, consisting of a set of droplets 
of triangular shape. In the valleys of the deposit, there 
are triplets of particles C with the structure shown in 
Fig. 3. Eventually, groups of four or more particles C 
may create such valleys, but they are much less probable 
then the triplets if p is small. These structures and the 
RSOS condition are responsible for the suppression of the 
growth process. 

For any p > 0, W s converges to a finite value with 
vanishing 1/L corrections. It contrasts to the behavior 
of moving phases, where the saturation width diverges 
as L a , with a > 0, in agreement with Eq.(2). Despite 
this remarkable difference on finite-size effects, W s will 
also be called saturation width here. Extrapolations to 
L — > oo give W s (p, oo) and the average saturation height 
H s (p, oo). The errors in H s (p, oo) are usually lower than 
1%, and the errors in W s (p, oo) are nearly 10%. 



II. NUMERICAL SIMULATIONS AND THE 
DYNAMICAL TRANSITION 

The main quantity of interest in deposition models is 
the interface width W of the deposit. In a surface of 
length L (L d columns), at time t, W is usually defined 



W(L,t) 



-1 1/2 



(1) 



where hi is the height of column i, the bar in h denotes a 
spatial average and the brackets denote a configurational 
average, i. e., an average over many realizations of the 
noise. 

In the pure RSOS model (p = 0), W obeys the dynamic 
scaling relation 



W re L a f (tL~ 



(2) 



The exponents a and z arepCpnsistent with the Kardar- 
Parisi-Zhang (KPZ) theorya, which provides a hydro- 
dynamic description of kinetic surface roughening. In 
d = 1, the KEjZ equation gives the exact values a = 1/2 
and z = 3/2E3. 

We simulated the RSOS version of the AC model for 
several values of p, most of them between p = 0.003 and 
p = 0.02. Substrates of lengths L from L = 256 to L — 
65536 were considered, with periodic boundaries. During 
the simulations, the time was measured as the number 
of deposition attempts per column, thus one time unit 
corresponds to L deposition attempts (accepted or not). 
For each p and L, we generated 10 sets of 10 3 different 
deposits each one, and calculated error bars from the 
fluctuations of the average values of the different sets. 

In all cases, the growth process fails at sufficiently long 
times, when the interface width attains a limiting value 
W s (p, L). In Fig. 2 we show a deposit for p = 0.1 and 
L = 128 in which no aggregation is possible. Notice 




FIG. 2. Example of a final static deposit for p — 0.1 and 
L = 128. 




FIG. 3. Triplet of particles C (filled squares), which oc- 
cupies most valleys of the static deposits, surrounded by A 
particles. 

In Fig. 4a we plot log H s (p, oo) x logp and in Fig. 4b 
we plot logW s {p, oo) x log p. Those quantities scale as 
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W s (p, oo ) ~ H s (p, 00) - p~ 



(3) 



with 77 = 1.509 obtained from the least squares fit of the 
H s data, and r\ = 1.515 obtained from the fit of the W s 
data. These relations show that the growth process will 
actually fail for any p > 0. 
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FIG. 4. (a) Saturation height H a (p, 00) in very large sub- 
strates versus probability p of incidence of particles C; (b) 
Saturation width Wsijp, 00) versus probability p. Solid lines 
are least squares fits. 



Our numerical results suggest the exact value 77 = 3/2, 
which can be obtained using scaling arguments, as fol- 
lows. The onset of triplets of particles C is responsible 
for the suppression of the growth process, and each block- 
ing configuration has probability of order p 3 . A mound 
of triangular shape (between valleys containing triplets 
of C) has height of order W s , then the number of parti- 
cles A in the mound is of order W s 2 . Thus, for small p, 



1/p 3 , giving 77 = 3/2. 



III. DYNAMIC SCALING 

The weak finite-size effects for large L suggest that 
a dynamic scaling relation in the static phase must be 
expressed only in terms of the probability p and the time 
t, while terms involving the length L will be (vanishing) 
corrections to scaling. 

For very large L, we propose the scaling relation 



w~p-y{t/r) 



p 



(4) 



where r is a characteristic time for the onset of corre- 
lations between the C triplets, and z is a dynamic ex- 
ponent, t is a measure of the number of layers of the 



deposit when these correlations appear, thus we expect 
that r ~ H s . Since H s also scales with exponent 77 (Eq. 
3), we obtain 



V : 



(5) 



In order to test relation (5) with the above exponents 
77 and z, we plot Wp v versus tp z in Fig. 5, considering 
three values of p: p — 0.005, p = 0.01 and p = 0.02. 
Those data were obtained in substrates with L — 65536, 
which are sufficiently large to minimize finite-size effects. 
The good data collapse in Fig. 5 confirms the validity of 
the scaling relation (5). 
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FIG. 5. Log- linear plot of Wp v versus tp z , with 
77 = 2 = 3/2, using data obtained in substrates with 
L — 65536 and probabilities p = 0.02 (squares), p = 0.01 
(triangles) and p — 0.005 (crosses). 



Finally, it is interesting to notice the divergence of the 
data for different p at t < 0.5p~ z , as shown in Fig. 5. 
At very short times, we expect that the interface width 
scales as in the pure RSOS model, with no dependence 
on p, because the effects of C particles are negligible. 
Then, the pure RSOS model regime, in which the width 
increases with time as i 1 / 3 , becomes just a transient re- 
gion for any p > 0. 



IV. DISCUSSION AND CONCLUSION 

We studied a deposition model with two types of par- 
ticles, A and C, in which incident particles can only stick 
at positions that have a neighboring A and if the RSOS 
condition is satisfied. For any flux of particles C, the 
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growth eventually fails, due to the RSOS condition and 
the formation of triplets of C. The saturation width W s 
is obtained in the static final configurations in sufficiently 
large substrates. Scaling arguments show that it scales 
as W s ~ p~ 3 / 2 for small p, and this result is confirmed by 
numerical simulations. The interface width W obeys a 
dynamic scaling relation involving the probability p and 
the deposition time t (Eq. 4). 

This model represents some growth mechanisms in the 
presence of impurities. As proposed in Refjj, it may de- 
scribe the effects of the deposition of an active particle B 
that reacts with a previously aggregated particle A and 
forms the inactive particle C. In the present RSOS ver- 
sion, small concentrations of the impurity may suppress 
the growth process, with the inactive particles forming 
the pinning centers. The blocking configurations depend 
on the particular model considered (for instance, they 
will change for different AH max), and the value of ex- 
ponent 77 depends on the number of particles C in those 
configurations. In a deposit with simple cubic lattice 
structure (which is more suitable for real applications) 
and AH max = 1, configurations with five particles C 
will form the pinning centers, and the supression of the 
growth process will also be observed. 

Previous works have also shown .transitions from a 
moving phase to a smooth phaseE3EJIi30. Usually, the 
roughening transitions are in the directed percolation 
(DP) universality class, but some of them possess other 
symmetries, e. g. the parity conserving class (PC). The 
smooth or anchored phases correspond to the active (or- 
dered) phases of DP, PC etc. In the moving phase, as 
the critical points are approached, the growth velocities 
continuously decrease to zero. The present model has 
many differences from those ones. First, the growth ve- 
locity changes discontinuously from a finite value at p = 
(pure RSOS model) to zero at p > 0. Furthermore,^ 
we consider the order parameters Mi defined in Ref.t-j 
(i = 1,2,.. .), we obtain Mj = in the static phase, since 
there is no preferential level for the pinning centers (see 
Fig. 2). Thus, this phase is not ordered in that sense. De- 
spite those differences, we expect that the analysis that 
led to the dynamic scaling relation (6) may be extended 
to other systems and may be useful to predict relations 
between growth exponents. 
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